Scaling behavior of linear polymers in disordered media 
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Folklore has, that the universal scaling properties of linear polymers in disordered media are well 
described by the statistics of self-avoiding walks (SAWs) on percolation clusters and their critical 
exponent i^sAw, with SAW implicitly referring to average SAW. Hitherto, static averaging has been 
commonly used, e.g. in numerical simulations, to determine what the average SAW is. We assert that 
only kinetic, rather than static, averaging can lead to asymptotic scaling behavior and corroborate 
our assertion by heuristic arguments and a renormalizable field theory. Moreover, we calculate to 
two-loop order i'saw, the exponent i/max for the longest SAW, and a new family of multifractal 
exponents u^'^K 

PACS numbers: 64.60.Ak, 61.25.Hq, 64.60.Fr, 05.50.+q, 05.70.Jk, 



In the past twenty years, the critical behavior of poly- 
mers in disordered media has generated a great deal of 
interest (for a recent review see The problem is 

relevant in a vast range of different fields. To name 
a prominent example, the transport properties of poly- 
meric chains in porous media might be exploitable com- 
mercially to enhance oil recovery. It has long been known, 
that polymers in disordered media are well modelled by 
self-avoiding walks (SAWs) on percolation clusters. The 
term SAW usually refers implicitly to average SAW. De- 
spite of many ideas put forward and extensive numerical 
efforts, the critical behavior of polymers in disordered 
media is still far from being completely understood. The 
most unsettling problems are, perhaps, on the analyti- 
cal side that the only existing field theoretic model for 
studying average SAWs, the Meir-Harris (MH) model [1], 
has trouble with renormalizability [1, 0] and that, on the 
numerical side, simulations lead to widespread results for 
the scaling exponent vsAW describing the mean length of 
average SAWs, see [H- 

A conceptual subtlety, that apparently has not been 
appreciated much hitherto, is the precise meaning of av- 
erage SAW. Viz. there are, essentially, two qualitatively 
different ways of averaging over all SAWs between two 
connected sites for a given random configuration of a di- 
luted lattice, one being static and the other being kinetic. 
In this letter, we conjecture that the statistics of linear 
polymers in disordered media has no asymptotic scaling 
limit, when static averaging is used. Since it is static 
averaging that has been commonly employed in numeri- 
cal work, many simulations might have suffered from this 
non-scaling behavior, which could explain the discrepan- 
cies between the numerical results for i'saw- In the fol- 
lowing, we first corroborate our conjecture by heuristic 
arguments. Then, we resort to renormalized field the- 
ory. It turns out that, for achieving renormalizability, 
one has to use kinetic averaging; static averaging leads to 



non-renormalizability. We then employ our field theory 
to calculate i^sAW, the corresponding scaling exponents 
for the shortest and longest SAWs and an entire family 
of multifractal exponents to two-loop order. Finally, we 
discuss the connection of kinetic averaging and the MH 
model. 

First, let us define what we mean by static and ki- 
netic averaging. In the framework of field theory, it is 
most convenient to model linear polymers in disordered 
media as SAWs between two connected sites x and y of 
a diluted lattice, were bonds are occupied with a prob- 
ability p, and to focus on the the length L{x,y) of a 
SAW (a random number proportional to the number of 
monomers of the corresponding polymer) rather than the 
Euklidian distance |x — y| of its endpoints 0. First, let 
us consider one given random configuration £ of the di- 
luted lattice. Averaging over all of SAWs belonging to 
the bundle B{x, y; £) of SAWs directed from y to x yields 
the mean length 



{L{x,y))e^K^\n 



(1) 



where L{'j) is the length of 7, p(7), with '^^p{"f) — 1, is 
a weight factor that depends on the averaging procedure, 
and K is the fugacity. Static averaging means that one 
simply uses ^(7) oc 1. Kinetic averaging, on the other 
hand, means that a SAW 7 earns a factor 1/z contribut- 
ing to 73(7) at each ramification where z — 1 other SAWs 
from the bundle B{x, y; €) split off. Experimentally rel- 
evant, however, is not {L{x,y))^ but rather its average 
[' ■ ']p O'^sr all configurations £ at fixed p subject to the 
constraint, that x and y are connected. This average is 
expected to exhibit scaling behavior, 

M(^,y) = [(L(x,y))4^|x-y|i/^--, (2) 

at a critical value Kc of the fugacity. 
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As we will demonstrate, SAWs on a percolation cluster 
are not merely standard fractals. Rather, they are mul- 
tifractals. In order to capture this multifractality, we de- 
fine the bond-weights uib = T,jeBix,y;C)Xb{l)ph) < 1, 
where Xbil) is one if the bond b belongs to the SAW 
7 and zero otherwise, and we introduce the multifractal 
moments 

L(")(x,2/)=^sfc™^ (3) 

b 

with Sb being the length of bond b. We will show that 
the scaling behavior of their quenched averages, 

M(")(x,2/)= [L(")(x,y)]^^|x-y|i/'''"', (4) 

is characterized by multifractal exponents v^"'' satisfying 
1/(0) = l/Dbb, v^^'' = vsAW, and iy^°°'> = i^, where Dbb 
is the fractal dimension of the backbone and v is the 
percolation correlation length exponent. 

It is well known Q that in a non-random medium 
(p — 1) the exponent t^sAW is the same for static and 
kinetic averaging. This may be not the case in a random 
medium, at least at the percolation point, and static aver- 
aging does not lead to a scaling law like Eq. ^ . Heuris- 
tically, this can be understood by employing the node- 
link-blob picture of percolation clusters in which a per- 
colation cluster connecting two terminal points, which is 
generically very inhomogeneous and asymmetric, can be 
envisaged as two nodes linked by tortuous ribbons that 
may contain blobs consisting of many short links in their 
interior. We will now use this picture to demonstrate 
that static averaging is unstable against coarse graining 
and that it therefore can not be expected to produce the 
correct asymptotic scaling behavior. Let us for simplic- 
ity consider the cluster sketched in Fig. [1] that features 
two links, one with and the other without a blob. With 
static averaging the (upper) link with the blob acquires 
a much larger weight then the other (lower) one even if 
it is much shorter than the link with the blob. Then, 
the statistics of the mean length is dominated by the 
short upper link with its many diS^erent SAWs induced by 
the blob. However, the weights change drastically upon 
coarse graining. Suppose we have some coarse graining 
procedure that culminates in condensing the blob into a 
single bond. After that, both links have the same weight. 
However, the lower one, since it is longer, now dominates 
the statistics. This demonstrates the instability of the 
weights of static averaging under real space renormaliza- 
tion as the group generated by repeated coarse graining. 
In contrast, kinetic averaging does assign the same weight 
to both links independent of the blob. Thus, kinetic av- 
eraging is stable under coarse graining even in a strongly 
inhomogeneous disordered medium. 

To fortify our arguments, we now turn to renormalized 
field theory. We will propose a theory for calculating 




FIG. 1: Percolation cluster in the node-link-blob picture. 

vsAW^ as well as the entire family z^("', that is renor- 
malizable, provided that kinetic averaging is used. This 
theory is based on the nonlinear random resistor network 
(uRRN) , where any bond on a d dimensional lattice is oc- 
cupied with a resistor with probability p or respectively 
empty with probability 1 — p. Our theory, is motivated 
by the well know fact, that the shortest and the longest 
SAW (the former is also known as the chemical path) can 
be extracted from the nRRN Jd| and its field theoretic 
formulation, the Harris model [7[ , by considering specific 
limits of the nonlinearity r of the generalized Ohm's law 
governing the bond resistors, 

~Vi^ p(^,j) l/.j-r sign/, J , (5) 

where Vi is the voltage at lattice site i, p(^ij^ is the 
resistance of bond (ij) and lij is the current flowing 
through that bond. As shown rigorously by Blumenfeld 
et al. [8], the shortest and the longest SAW correspond 
to r — > -1-0 and r — > —0 respectively. Evidently, M(a;, y) 
must lie in between the average length of the shortest 
and the longest SAW, which are, of course, very differ- 
ent. Since M{x,y) sits somewhere in this discontinuity 
at r = 0, it is not known how to extract it from the 
nRRN by a limit taking. Therefore, we propose here to 
study the average SAW by using our real-world interpre- 
tation [9, 10, 11, 12], in which the Feynman diagrams for 
the nRRN are viewed as being resistor networks them- 
selves. The idea is to put SAWs on these diagrams. That 
this idea is fruitful can be checked explicitly at the in- 
stance of the chemical path. Our approach reproduces 
to two-loop order the corresponding exponent Vm\n well 
establish from dynamical percolation theory [l^. 

Our field theory is based on the Harris model as de- 
scribed by the Hamiltonian 

u = j d'^x5]{i^2+i(v^)2-K|^(-9;)^+v+fv'^}, 

e 

(6) 

where is a replicated discretized voltage taking 
on {2N)^ values on a ZJ-dimensional torus: 6 = 
jj{ni,--- ,nD) withrij = ~N + l,-N + 2,- ■ ■ ,N-l,N. 
Lp — y'(x, 6) is the order parameter field, a continuum 
analog of a Potts spin. It transforms according to the 
one irreducible representation of the symmetric (permu- 
tation) group S(2N)'^ and thus, the model features only 
a single coupling constant, g. r and w are strongly rel- 
evant critical control parameters. The scaling behavior 
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FIG. 2; (Color online) Static and kinetic rule for averaging 
SAWs. 



of SAWs is associated with the renormalization of w in 
the replica limit £) — > 0. For details on the Harris model, 
we refer to 0, 0, The diagrammatic perturbation 

theory of the Harris model can be formulated in such 
a way that the Feynman diagrams resemble real RRNs. 
In this approach, which we refer to as real world inter- 
pretation, the diagrams feature conducting propagators 
corresponding to occupied, conducting bonds and insu- 
lating propagators corresponding to open bonds. The 
conducting bonds carry replica currents conjugate to the 
replica voltages 9. The resistance of a conducting bond is 
given by its Schwinger parameter In the following, 
we will refer only to these very basic aspects of the real 
world interpretation which will be sufficient to follow the 
main line of argument. For further details on the real 
world interpretation, see Refs. [l^, 11. 

Now we employ the real world interpretation to study 
the average SAW. As an example, let us consider the two- 
loop diagram that resembles the node-link-blob cluster 
in Fig. [T] and where all internal propagators are conduct- 
ing. For determining the contribution of this diagram 
to M{x,y), the essential step in our approach is to find 
out the length L of an average SAW on that diagram. 
We can either apply the static or the kinetic rule to cal- 
culate L according to Eq. Q with a = 1, where now 
bond b is replaced by conducting propagator i and the 
Schwinger parameter 1^ Si is interpreted as the corre- 
sponding length. As visualized in Fig. (I2|), static and 
kinetic averaging yields different results. The static rule 
gives L^** 



(si -|- S4 -I- S5)/3 -|- 2(s2 + S3)/3, whereas the 
kinetic rule produces L^^'") = (si-|-S2-l-S3)/2-|-(s4-|-S5)/4. 
The remaining steps in calculating the diagram essen- 
tially textbook matter IJ]. It turns out, that the static 
rule does not lead to a renormalizable theory. The rea- 
son is easily shown. In the 2-loop calculation of our 
example diagram, the non-primitive divergencies aris- 
ing from the sub-integrations of the 1-loop self-energy 
insertion must be cancelled through the counter-terms 
introduced by the renormalization of this 1-loop inser- 
tion. However, the weights of L'^*'*) are not in confor- 
mity with the weights arising in the corresponding 1- 
loop diagram with the counterterm insertion: crunching 
the insertion to a point (corresponding to S4 + S5 0) 
leads to i^''*) = si/3 + 2(s2 + S3)/3 in contrast to 
2^ (kin) ^si/2 + (s2 + S3)/2, which is equal to L of the 1- 
loop self-energy diagram with a point insertion. Hence, 
only the weighting according to the kinetic rule works 
correctly in that it leads to a cancellation of non-primitive 
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FIG. 3: (Color online) The e-expansions of the exponents Vmin 
(blue upper curve), i^sAw (red middle curve), and i/max (green 
lower curve). Possible extrapolations at low dimensions d are 
shown by broken lines. The squares with error bars symbolize 
numerical results for us asm as compiled in [ll]. 



divergencies by one- loop counterterms. Thus, the static 
rule has to be rejected on grounds of renormalizability. 

Besides revealing the imperative of kinetic averaging, 
this theory yields two-loop results for the SAW exponents 
vsAW and t'max, which previously have been calculated 
(correctly) only to one- loop order 0, 0| , and the family 
which is entirely new: 
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where e — 6 — d. vsAW is given by vsAW = i^''^''- Our 
result for z^sAW is compared to the available numerical es- 
timates, to our result for the longest SAW, and to the well 
known exponent z^min for the shortest SAW f9l, [13, [l3| in 
Fig. The following points are worth noting: (i) v^"^ 
does not depend on a in a linear or affine fashion which 
implies that SAWs on percolation clusters are mulitfrac- 
tal. (ii) u^"^ is in absolute agreement with the well known 
results for Dbb and ly in the cases a = and 00, respec- 
tively, (iii) Vmin and z^max are not related to the family 
i/^"-'. (iv) the theory is renormalizable for arbitrary a if 
and only if kinetic averaging is used. 

As mentioned above, the usual framework to study av- 
erage SAWs on percolation clusters is the MH model as 
described by the Hamiltonian 



(9) 



Here, ^-fe = {^fe;al W}, 1 < fc < n is an or- 
der parameter field conjugate to a n-fold replicated m- 
component Heisenberg spin with vector-indices ii running 
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from 1 to m and replica indices a/ G {1, . . . ,n} ordered 
so that ai < • ■ • < ak- The n repUcas transform ac- 
cording to n different irreducible representations of the 
direct product of the symmetric group S'„ and the or- 
thogonal (rotation) group SO{m). — 'Y^^wik} , and 
is a symbolic notation for the sum of the products 
of three '^k fields. Only those cubic terms are allowed 
for which all pairs {i,Oi) appear exactly twice. In this 
model, one can extract j^saw from the renormalization 
of the relevant control parameter wi upon taking the 
replica limit n — > 0. The MH Hamiltonian ^ is non- 
renormalizable as it stands. One difficulty, that has been 
pointed out by Le Doussal and Machta I3| several years 
ago, is that the critical values {r^} of the control param- 
eters are different for different fc, i.e., the model is highly 
multicritical. A second problem, that to our knowledge 
has not been discussed hitherto, is that the order pa- 
rameter fields ^fe for different k belong to different ir- 
reducible tensor-representations of underlying symmetry 
group, Sn X SO{m), see above. Hence, strictly speak- 
ing, one needs independent coupling constants gk,k' ,k" 
for each product '^k ^fe' ^fe" [note that this is (i) not im- 
plemented in the original MH Hamiltonian ([9]) and (ii) 
different in the Harris model], and the fields need k- 
dependent renormalization factors [l5l |. Recently, these 
difficulties caused the failure of a two-loop calculation of 
vsAW by von Ferber et al. Q. 

As far as its application to the average SAW is con- 
cerned, the renormalizability of the MH model can be 
rescued by a specific interpretation of the replica limit 
which has close ties to kinetic averaging. Our analysis of 
the MH model (details will be given elsewhere [17]) led 
to the following key findings: If the replica limit is taken 
after all summations over all possible arrangements of in- 
ternal replica indices of a diagram, then the MH model 
reproduces static averaging and, as demonstrated above, 
is not renormalizable. If, however, the replica limit is 
taken, in the spirit of Ref. [3], as early as possible, i.e., 
loop after loop, or at least for each renormalization part, 
then the MH model reproduces kinetic averaging. This is 
the only interpretation of replica limit of the MH-model 
that leads to a renormalizable theory of SAWs in disor- 
dered media. With this interpretation, the MH model, 
in particular, produces the same result for i^sAW as our 
real-world interpretation and thereby provides an impor- 
tant consistency check for the validity of the application 
of the latter to SAWs. 

Closing, we would like to emphasize that our renor- 
malization group arguments are, although certainly well 
founded, not rigorous in the sense of a mathematical 
proof since they rely on our real world interpretation of 
Feynman diagrams. This interpretation thrives on anal- 
ogy and there exist to date no rigorous mathematical ar- 



guments on how far its validity extends. However, given 
all its successes in the past, we would be surprised if it 
failed in describing SAWs on percolation clusters. The 
well known MH model, when interpreted carefully, cor- 
roborates the imperative of kinetic averaging and con- 
firms our two-loop result for j^saw- 
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